Quantum jumps induced by the center-of-mass motion of a trapped atom 
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We theoretically study the occurrence of quantum jumps in the resonance fluorescence of a trapped 
atom. Here, the atom is laser cooled in a configuration of level such that the occurrence of a quantum 
jump is associated to a change of the vibrational center-of-mass motion by one phonon. The statistics 
of the occurrence of the dark fluorescence period is studied as a function of the physical parameters 
and the corresponding features in the spectrum of resonance fluorescence are identifled. We discuss 
the information which can be extracted on the atomic motion from the observation of a quantum 
jump in the considered setup. 

PACS numbers: 



I. INTRODUCTION 

Quantum jumps are one of the earliest concepts in 
quantum mechanics, which were introduced to denote the 
atomic transitions between electron bound states due to 
photon absorption or emission as predicted by Bohr [l|. 
This terminology is now used to designate the observa- 
tion of a sudden change in the fluorescence signal of laser- 
driven atoms, which is associated with the occurrence of 
transitions between electronic states 0-0] • Since their 
first observation [2tJ], quantum jumps have been exten- 
sively discussed in the literature. For an overview we 
refer to the reviews [5|-|7[. Quantum jumps have inspired 
theoretical descriptions and methods in quantum optics, 
such as quantum trajectories [g, I8l-[l0|. and are at the 
basis of several theoretical proposals for post-selected 
quantum state preparation, some examples are found 
in Ref. |lll - [l3| . and for detecting the quantum state of 
trapped atoms [13, [3| • Most recent experiments used 
quantum jumps in order to demonstrate single photon 
absorption by single trapped atoms [15| , and analyzed 
their statistics in order to demonstrate quantum corre- 
lations between the absorbed photon and a second, en- 
tangled photon revealed at a detector [la|. Further ex- 
perimental work used quantum jumps in the signal mon- 
itoring the number of photons inside a resonator as the 
basic element for inrplementing feedback loops, restoring 
the quantum state of a cavity mode [l7[ . These milestone 
experiments demonstrate once more that quantum jumps 
can be a powerful tool for monitoring and manipulating 
the quantum state of physical systems. 

Quantum jumps have been extensively discussed in 
connection to the preparation and measuring of the quan- 
tum state of a trapped atom p, [Il|, [ij, [I^ . In a suit- 
able setup their analysis can provide information on the 
dynamics and the state of the atomic center of mass de- 
gree of freedom. Information on the quantum state of 
the atoms can be also extracted by the spectrum of reso- 



nance fluorescence, which, when taken at the steady state 
of laser cooling, reveals the relevant scattering processes 
characterizing the dynamical steady state [19|, |20| and 
can provide a measurement of the coo ling rate and steady 
state phonon number occupation |2l| - [2a |. 

For the detection of quantum jumps one usually mon- 
itors the fluorescence light of a strongly driven atomic 
transition. When due to a additional weak coupling the 
electron jumps to a third, metastable state, the sud- 
den stop of the fluorescence reveals the occurrence of 
the jump. The properties in the resulting spectrum of 
resonance fluorescence have been theoretically analyzed 
in Ref. ^26.] and measured in [23]. General features in 
the emission and absorption spectrum of three-level sys - 
tems have been theoretically investigated in Refs. [28i,i29|, 
and experimentally confirmed j30l . i31i] . During the last 
decades, three-level systems have become a paradigmatic 
system for observing quantum interference |32l434| . Sig- 
natures thereof can manifest themselves in the spectrunr 
of resonance fluorescence [2^, I35l-l37ll . 



In this work we address the question what are the spec- 
tral properties associated with quantum jumps, carrying 
information on the center-of-mass motion of a trapped 
atom. In particular, we establish the connection be- 
tween the time dependent fluorescence and the corre- 
sponding spectrum originating from the quantum mo- 
tion of the atomic center of mass. For this purpose, we 
analyze a setup, constituted by a single trapped atom 
driven by laser fields which cool its center-of-mass mo- 
tion, such that the occurrence of quantum jumps in the 
atom's fluorescence is associated with a change of its 
center-of-mass motion by a single vibrational excitation. 
The configuration of electronic levels and the setup are 
sketched in Fig. [1] The telegraphic noise, detected in 
the fluorescence, is analyzed both in time and in fre- 
quency, thereby providing an extensive description of the 
quantum jump process induced by the quantum fluctu- 
ations of the atom's center-of-mass motion. We identify 



and characterize the spectral features associated with the 
quantuni jumps, which reveal a change of the atomic mo- 
tional state by one vibrational excitation. Our analysis 
shows that the spectral components connected with the 
motion-induced quantum jumps are proportional to the 
ratio between the dark and bright time scales of the in- 
termittent fluorescence, and scale with the temperature 
of the laser-cooled atom. 

The article is organized as follows. In Sec. [TT] the setup 
and the corresponding theoretical model are introduced. 
In Sec. mil the time scales of bright and dark periods in 
the resonance fluorescence are evaluated using the de- 
lay function [3, [2^|, and the basic steps leading to the 
evaluation of the spectrum of resonance fluorescence, are 
summarized. Moreover, the results are discussed in this 
section. In Sec. HV] the conclusions and a discussion on 
possible outlooks are presented. The appendices report 
the details of the calculations at the basis of the results. 



II. THE SYSTEM 

We consider an atom of mass M whose center-of-mass 
motion is confined by a harmonic trap. The atom is 
driven by lasers, and the relevant electronic levels, which 
are involved in the dynamics, are the stable state \g) and 
two excited states |1), |2) at frequencies ui and uj2, re- 
spectively. The level scheme is sketched in Fig. [T] a). 
States |1) and |2) are excited by means of lasers and de- 
cay radiatively into the ground state \g) with rate 71 and 
72, respectively, where we assume that 72 ^ 71. This 
configuration has been considered for observing quantum 
jumps in atomic systems: a long-living excitation of state 
1 2) leads to interruption of the fluorescence along the 
transition \g) o |1) [2^|33- In the present setup, how- 
ever, there is no excitation along the transition \g) <->■ |2) 
if the atom is a point like particle, since the laser field 
driving this transition is a standing wave and the trap 
center coincides with one of its node, see Fig.[T]b). Hence, 
only the mechanical effects of light will induce dark pe- 
riods in the resonance fluorescence signal. 

We now discuss the mechanical effects of the light on 
the atomic motion. In the considered scheme, the laser 
driving the transition \g) ^ |1) is a running wave with 
frequency wli and wavevector fci, which Doppler cools 
the atomic motion [40J. The laser driving the transition 
\g) O |2), instead, is a standing wave with frequency wl2 
wave vector ^2, where one of its nodes coincides with 
the center of the trap, as shown in Fig. [T]b). Excita- 
tions of the state |2) hence occur due to fluctuations of 
the atom's position about the center of the trap. In this 
work we consider that the atomic motion is in the Lamb- 
Dicke regime, such that the extension of the center-of- 
mass wave packet is much smaller than the laser wave 
length and the mechanical effects are hence a small cor- 
rection to the motion. In this regime, interruptions in 
the fluorescence on the transition \g) o |1), quantum 
jumps, will then occur due to a mechanical excitation or 
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FIG. 1: A single atom is confined in a harmonic trap and 
laser-cooled, a) Relevant levels, showing the dipole transi- 
tions \g) -f-)- |1) and \g) -f-)- |2) which are driven by lasers with 
Rabi frequency Jli and ^2, respectively. The state \g) is sta- 
ble, and excited states |1), |2) decay into it with rates 71 and 
72. In this paper 72 ^ 71. b) Setup: The harmonic po- 
tential with frequency z/, confines the atomic center of mass 
motion which is assumed to be along the a;-direction. The 
laser couphng the transition \g) -f-^- |1) is a running wave at 
wave vector ki forming the angle (j>i with the 2;-axis, as illus- 
trated by the black arrow in the figure. The laser couphng 
the transition \g) ^-s- |2) is a standing wave at wave vector 
^2 forming the angle (/!>2 with the x- axis. In the figure it 
is represented by the shaded region with sinusoidal envelope. 
The trap center coincides with one node of the standing wave. 
A detector collects the photons emitted at an angle tp with 
the axis. Inset: Temporal behavior of the resonance fluores- 
cence, consisting of bright and dark periods at time scales 
Tb and To , respectively. A quantum jump to a dark period 
corresponds to a transition to state |2), and is due to atomic 
fluctuations around the trap center at a node of the intensity 
of the laser coupling to state |2). 



de-excitation of the atomic center of mass motion. 



A. Theoretical model 

The dynamics of the density operator g, describing the 
internal and external degrees of freedom of the atom, is 
governed by the master equation 



dg 



£g, 



(1) 



where the Liouville operator C is defined as 



1 



Cq = —: \H, g] + JCg, 
in 



with the Hamiltonian H, giving the coherent evolution of 
the system, and the relaxation superoperator /C describ- 
ing spontaneous emission. In detail, the Hamiltonian can 
be decomposed into the sum of the terms 



H^H, 



Hu 



V, 



consisting of the center-of-mass motion in the trap, which 
we assume here to be along the x— axis, 



i^cxt ^hiy{a'' a +1/2), 



with the annihilation and creation operators a and a^ of 
a vibrational excitation at energy hi'. Furthermore, it 
contains the electronic states' energy, 

i:?i„t-Mi|l)(l| + M2|2)(2|, 



together with the sufficient condition rjj^J {n) <^ 1, where 
(n) = {a^ a) is the mean number of phonons. The mean 

(2) number of phonons (n) is reached by an interplay of the 
mechanical effects of both lasers. The laser driving the 
transition \g) <->■ |1) is tuned to the red side of the atomic 
resonance, 6i = 71 /2, and performs Doppler cooling on 
the atom. In absence of the second laser, the photons 
scattered along this transition, contain the information 
about the steady state of Doppler cooling. We also as- 

(3) sume that the second laser realizes sideband cooling, such 
that absorption of a photon along the transition \g) -^ |2) 
corresponds with largest probability to a de-excitation of 
the atomic motion by one phonon |4C)I443 | . Although it is 
not essential for our study, we have optimized the param- 
eters in order to minimize the final temperature of the 
combined laser cooling. A brief discussion of the cooling 
dynamics and of the choice of parameters for optimizing 
the steady-state temperature is reported in App. |B] Be- 
low we introduce the expansion of the Liouvillian C in 

(5) powers of the Lamb-Dicke parameters. 



(4) 



which are reported in the frame rotating at the frequen- 
cies of the lasers, such that 5j = ujj—uji^j. The laser-atom 
interaction in electric dipole approximation has the form 



V =h— [|l)(g|e*'=i^™'''*i 

1^2 



H.c] 



+h^ sin (k2X cos (j)2) [|2) (<?| + H.c] , (6) 

with the position x — ^{a -\- a^) of the atomic center of 



mass motion, and ^ = \ ^m^- The interaction includes 
the mechanical effects associated with the absorption and 
emission of a photon of the laser, which are weighted 
by the projection kj cos (j)j of the lasers' wave vectors on 
the X-axis. Finally, the superoperator describing sponta- 
neous decay takes the form 

^^ = E f (2|5>01e,b-)(<?l - q\j){3\ \3){3\q) , (7) 
i=i,2 



B. Lamb-Dicke expansion of the Liouvillian and 
spectral decomposition 

Master equation ([2]) can be formally solved using the 
spectral decomposition of the Liouvillian £, namely, the 
left- and right eigenelements at eigenvalue A obeying 



Ct = Xg^, g^C = A^ 



^A 



(10) 



Orthonormality and completeness are defined with re- 
spect to a scalar product given by the trace, such that 
Ti{g'^g'^ } — (5a, A'- Correspondingly, the projectors V^ 
can be introduced, whose action on an operator X is 



r^X = g^Tr{g^X} 



(11) 



such that the completeness relation can be expressed as 



j2p. 



1. 



(12) 



where 



gj = / auw[u)e ' ge ^ . 

-1 



(8) 



The weight function w{u) is the normalized dipole radia- 
tion pattern, giving the probability of emitting a photon 
at an angle u = cos 9 with the motional axis. In the fol- 
lowing we will take 'w{u) = 3(1 -|-u^)/8, corresponding to 
the transitions \g) <->■ |1),|2) being characterized by the 
selection rule for the angular momentum A J^ = ±1. Fur- 
thermore, we will assume the Lamb-Dicke regime, corre- 
sponding to taking the Lamb-Dicke parameter rjj — ^kj 
such that it fulfills the relation |41j] 



Vj = 




«1, 



(9) 



We note that the steady state is a right eigenelement of 
the Liouvillian C at eigenvector A = 0, while the corre- 
sponding left eigenelement is g^^'^ — 1 [43[. 

Eigenvalues and eigenvectors can be evaluated using 
the Lamb-Dicke expansion. The Lamb-Dicke expansion 
of Liouvilhan C, Eq. ©, up to second order gives £ — 
Cq+Ci +C2 -l- . . . , where the subscripts label the order in 
the expansion and the individual Liouvillian terms take 
the form [Hi 

-= ^[^cxt, q] + (jf^iHint + Vo, g] + ICog] , (13) 



Cig^ -r[Vix,g], 
in 

C2g^ —-[V2x'^,g] +IC2g, 
2in 



(14) 
(15) 



for the zero, first and second order terms, where we have 
used 



K = 



d'^Vix) 



9a;" 



(16) 



x=0 



and 



K,oQ= Y. y(2|.9)0>b-)(<?|-e|j>01-|j)0» 



J = 1.2 



l]k 



1 

where W2 — J du'w{u)u'^. 
-1 
In zero order, the internal and external degree of free- 
doms are decoupled, as one can also see in Eq. ([T3|) : Liou- 
villian Ce describes the equation of motion of a harmonic 
oscillator, while Liouvillian Ci gives the internal dynam- 
ics of an atom at rest. Correspondingly, in zero-order 
perturbation theory the eigenvalues read 



Ao — Ae + Ai 



(17) 



where Ae = i)^i^ {i = 0, ±1, ±2, . . . ) are the eigenvalues 
of £e and Ai are the (nine) eigenvalues of £1. The steady 
state, right eigenvector at Aq = 0, takes the form 



9f 



Pstl^st 



(18) 



with the density matrix pst defined for the internal de- 
grees of freedom and obeying Cipst = 0, while /igt is 
the density matrix for the center-of-mass degrees of free- 
doni, which solves the effective equation for laser cool- 
ing [22, [231 ■ Projectors on the corresponding subspaces. 



K 



/ E ' I ' 



(19) 



are a product of externals projectors 

V^^^'^''X = Y\n){n\X\n + l){n + l\ (20) 

n 

and internal projectors V^ ' of £1. The internal projectors 
Vi are constructed by the left- and right eigenelements 
to Ai, Pi^ and /5^^, of £1, see App. [X] Note that the 
external subspaces are infinitely degenerated. Higher or- 
der projectors Vi, V2,- ■ ■ and eigenvalues Ai, A2,. . . can 
be calculated using standard techniques of perturbation 
theory and can be found in [1^, [23] . It turns out that 
Ai = and VqCiVq = since £1 does not couple states 
of subspaces with same eigenvalue Aq [1^; ^M ■ 



III. QUANTUM JUMPS AND SPECTRUM OF 
RESONANCE FLUORESCENCE 



In this section we identify the characteristic time 
scales, determining the duration of the bright and dark 
periods in the resonance fluorescence as a function of 
time, and then determine the spectrum of intensity of 
the emitted light. 



A. Quantum jumps: Dark and bright periods 

We now focus on the fluorescence signal as a function 
of time. We show that, studying light scattering by the 
atoms at steady state of laser cooling, the fluorescence 
exhibits bright periods, due to light scattering along the 
transition \g) ^ |1), and dark periods which are due to 
the shelving of the electron in the state |2). The transi- 
tions to the state |2) can only take place due to the finite 
extension of the motional wave packet of the atom at the 
node of the laser standing wave. The corresponding pro- 
cesses are therefore of higher order in the Lamb-Dicke 
expansion and associated with a change of the atom's 
center-of-mass state. In the Lamb-Dicke regime, the time 
scales characterizing the bright and dark periods of the 
fluorescence are well separated: Denoting by Tb and Tq 
the time scales of bright and dark periods, respectively, 
thenrs>TD. 

In order to evaluate Tb and Tq we consider the 
waiting-time distribution P{t), which is the probability 
that no photon emission occurs in the interval of time i, 
provided that a photon was detected at the instant i = 0. 
The waiting-time distribution is defined as [y, |8| 



P{t) ^ Tr{e-'^'=«*/'"Vffffe'-f^°«*/'^}, 



where 



ffni 



H 






(21) 



(22) 



is a non-Hermitian operator, with H given by Eq. ([S]). 
and pgg is the density matrix describing the state of the 
system, assuming that at t — the system is found in the 
internal ground state. In our case the density operator 
Pgg corresponds to performing a measurement at time 
t = on the steady state density matrix, whose outcome 
is the projection onto state \g). Hence, in our formalism 

Pgg = l.9)(5|TriPst = \9}{9\Pst , 

where /igt is the steady state density matrix, resulting 
from the laser cooling process (see App. |B|). Figure [H 
displays the waiting function P{t) as a function of time 
for a specific choice of the parameters. One observes 
a fast decay of P{t) for short times, which is followed 
by a slow decay with a smaller rate at larger times. A 
time scale t can be identified, which separates the two 
behaviors. From P{t) we can estimate the average time 
scales 

Tb = [1 - P{t)]-' |t + P{t)-^ J^^ dt'P(t')^ , (23) 

dt'P{t') (24) 

for the bright and dark periods, respectively. We cal- 
culate the two time scales perturbatively up to second 
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FIG. 2: Waiting-time distribution P{t) as a function of the 
time interval between two subsequent quantum jumps events, 
in units of 1/72- The waiting-time distribution corresponding 
to the solid line has been evaluated numerically from Eq. (|21|) . 
whereas the dashed curve follows from the perturbative treat- 
ment in App. [C] Two well separate time scales can be ob- 
served, from the short one follows the bright period Tb, the 
long one corresponds to the time scale of the dark period, 
Td- For the considered parameters, they are in agreement 
with the expressions given in Eqs. (|25p and (|26l) . The pa- 
rameters are 71 = 12;^, 5i = 6v, Qi = 2.5z^ and 72 = 0.015;/, 
^2 ~ 0.5i/, 52 = 0.87^. The Lamb-Dicke parameters are taken 
to be rii = ri2 — 0.05, the angle of the laser wave vectors are 
01 = it/9 and (f>2 = 0, the detector is at angle tp = 47r/5 (not 
relevant for this plot). The mean phonon number at steady 
state is (n) « 0.3. 



order in the Lamb-Dicke parameter. The details of the 
calculation are reported in App. [Gland yield 



Tr 



Tn 



1 



AT' 



7y| cos^ 
1 

72' 



m 



j.iO) 



where 



T(°) = (71P11)" 



(25) 
(26) 

(27) 



is the time scale of photon emission from state 1 1) , being 
pii the occupation of state |1) at steady state and leading 
order in the Lamb-Dicke expansion, pn = ■^s/{s + 1), 
with the saturation parameter [3^ 



nl/2 



(28) 



for the transition \g) o |1). The explicit form of 
the temperature-dependent parameter F is given in 
Eq. ()C12|) . We assume small saturation s ^ 1 and 
72 ^ S71/2, corresponding to the assumption that the 
rate 7^ = s7i/2 of incoherent scattering on the transi- 
tion \g) <-> |1) is faster than the decay rate of state |2). 
In this limit, we find 



7i 
4(n) 



(29) 



Using Eq. ([29]) and Eq. ([271), for s < 1, in Eq. ([25 
scale for the bright periods takes on the form 



Tb 



-ql cos^ 



h{n] 



7( 



the 



(30) 



which corresponds to the time scale for scattering a pho- 
ton along the transition \g) <->■ |2). Moreover, since scat- 
tering along this transition is exclusively due to the me- 
chanical effects of light, the time scale of the bright pe- 
riods becomes larger, the smaller are either the atomic 
spatial fluctuations about the trap center, since they scale 
with Ax^ ~ ^^ cos^ (t>2 {n) , or the smaller is the mechani- 
cal effect of the standing wave laser on the atom (which 
indeed vanishes when the laser propagates perpendicu- 
larly to the atomic motion, cos 02 =0)- In particular, 
when cos 02 = there are no excitations of state |2) and 
hence no quantum jumps are observed, corresponding to 
Tb — >• 00. 



B. Evaluation of the spectrum of resonance 
fluorescence 



The relation between the temporal behavior of the fluo- 
rescence signal, governed by the dark and bright periods, 
with the spectrum of resonance fluorescence was studied 
in Refs. [y, [2y] for an atom at rest. There, the authors 
found the occurrence of a sharp peak at the laser fre- 
quency, whose width scales with the inverse lifetime of 
the shelving period. We now extend these studies to the 
case of a trapped atom in the Lamb-Dicke regime, using 
the formalism first introduced in [2lL i22i] and then fur- 
ther developed in [13, 120, 1231 ■ We evaluate the spectrum 
of resonance fluorescence by decomposing the spectrum 
in the eigenbasis of the Liouville operator, Eq. ^ and 
evaluate the resulting expressions perturbatively in the 
Lamb-Dicke parameter. Such an approach allows us for 
discussing the relevant components of the spectrum be- 
longing to the two transitions \g) <->■ \j) individually. 

We consider the spectrum of the fluorescence light scat- 
tered by the atom along the transition \g) o |j) and 
recorded in the far-fleld. The atom is assumed to be in 
the (dynamical) steady state of laser cooling. The com- 
ponent of the spectrum at frequency uj is proportional to 
the function 

S'^'\oo) = Re f dte~'^^-^^^^'{D<^^^\t)D(J'>),t (31) 



with 



D — \9){j\exp[—ikjXC0S->p], 



(32) 



giving the positive- frequency component of the dipole op- 
erator, which describes atomic de-excitation and the re- 
coil of the emitted photon which is measured by the de- 
tector [4^ at angle ^p. Note that (•)st = Tr{-gst}, where 



Qst is the steady state defined by Cgst = 0. In Eq. ((3T|) 
we omit an overall prefactor including the dipole radia- 
tion pattern of the transition. In the following discussion 
we omit the superscript j when discussing the general 
properties. 

Using the quantum regression theorem, we can rewrite 
the expectation value in Eq. (PT|) as {D^{t)D)st — 
Tr{D^e'-'^Dgst}- Inserting the completeness relation P^ 
between the operators, the power spectrum, 
takes the form 



Eq. dH 



5H=Re^. 



1 



i(a; — wl) ~ A 



F(A) 



where 



F(A) = Tr{D^V^Dg,t} 



(33) 



(34) 



is a weight function of each spectral contribution. The 
explicit form of the spectral component can be deter- 
mined using the Lamb-Dicke expansion of the Liouville 
operator and of the Dipole operator D{x). In particular, 
in zero order in the Lamb-Dicke parameter the spectrum 
has the form 



5o(a;) = Re^ L__Fo(Ai), 

^ i{ijj - wl) - Ai 



(35) 



with i^o(Ai) = TrjDj-Poi'oPst} and Do = D{{)). Equa- 
tion (l35t gives the spectrum of the bare atom at rest. It 
consists of the 5-like elastic peak at the laser frequency, 
corresponding to the eigenvalue Aj = 0, and of the inelas- 
tic components, corresponding to Ai 7^ 0. For a two-level 
transition it corresponds to the spectrum evaluated by 
Mollow [13. 

The first non- vanishing contributions to the spectrum 
due to the mechanical effects of light are at second order 
in the Lamb-Dicke expansion, and we denote them by 
the term 



52(w) = Re^ 



1 



i{uj — LUl) — \ 



F2{X), 



with 



F2(A)= Y. Tr{i?tp^A^^^f}, 

a+[3+-,+S=2 



(36) 



(37) 



and Dn ~ d'^D{x)/ dx'^\x=o/n\ is the corresponding coef- 
ficient of the expansion of the dipole operator. For later 
convenience we split the second order spectrum S2{ijj) 
into the correction to the elastic peak, corresponding to 
the contributions in second order in the Lamb-Dicke pa- 
rameter at frequency lo — ui^ and here denoted by 5*01 (w), 
the contribution at the frequencies w = wl ± i^, corre- 
sponding to the motional sidebands of the elastic peak, 
given by the term 



5sb(t^)=Re Y. 



1 



i(cj — cjl) — A 



^2(A), (38) 



and finally the contribution from inelastic scattering, 

1 



5'inGi(w) = Re y^ 



i(aj — wl) — A 



^2 (A), (39) 



such that S2(oj) — Sc\(io) + S'sb(w) -I- S'inci(w). In general, 
the inelastic component at second-order in the Lamb- 
Dicke parameter is a small correction to the inelastic 
component at zero order. Detailed examples on how to 
evaluate the spectra can be found in [1^ [2j] . 

The spectra are solved evaluating numerically the 
traces Fq and ¥2- Nevertheless, in some limit analyti- 
cal expressions can be derived for the terms of interest 
to this work, namely, the ones related to the appearance 
to quantum jumps in the temporal behavior of the reso- 
nance fluorescence. 



C. Spectra of resonance fluorescence: Results 

Figure [3] a) displays the spectrum of resonance fluo- 
rescence of the light emitted from the atom along the 
transition \g) -f-)- |1), where the elastic peak, at w = wli, 
is not shown. Three peaks are visible: (i) the two side- 
bands of the elastic peak, corresponding to the photons 
emitted at frequency wli -t- i^ and wli — i^, and which 
are associated with the scattering processes which cool 
and heat, respectively, the center-of-mass motion by one 
phonon, and (ii) a central peak at w = wli, which is 
due to the inelastically scattered photons. These three 
features are in second order in the Lamb-Dicke parame- 
ters, and are hence solely due to the mechanical effects of 
light. Further components, like for instance the inelastic 
components at zero order in the mechanical effects, are 
nonvisible and completely overtopped by the signals due 
to the mechanical effects. The sidebands of the elastic 
peak have the form already known from Doppler cooled 
atoms in the asymptotic stage of cooliiig, and have been 
object of several extensive studies |19l. [2l| - |23{ . Here, we 
focus on the central peak, which is a novel feature with 
respect to the spectrum recorded by a two-level atom 
which undergoes Doppler cooling. An analytical form of 
the behavior of the spectrum around this point can be 
extracted from our theory. For s ^ 1 and assuming 72 
to be the smallest decay rate, (72 ^ ■S71/2) it has the 
form 



s^A^) 



72 



Td s 



(w- WLi)^ -}-7| Tb 2' 



(40) 



Al = 0,AE = ±i 



where the details of the calculations are reported in Ap- 
pendix |DJ The analytical result in a more precise form, 
see Eqs. (|Dip and (jDSp . is compared to the numerical 
one in Fig.[3]b) and shows good agreement. This contri- 
bution is the one giving rise to the pedestal of the elastic 
peak due to the inte rrup tion of fluorescence, and was 
first described in Ref. [2g| for the case of an atom at rest. 
This central peak can be understood as a consequence 
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FIG. 3: a) Spectrum of resonance fluorescence of the light 
emitted by state |1). The inset shows the magnified sidebands 
of the elastic peak. The visible peaks are due to mechanical 
effects of light, b) Magnification of the central narrow peak. 
The solid line represents the spectrum from Eqs. (|35[1 and 
P6|l . while the dashed line corresponds to the approximation 
using the expressions (|D1|) and (|D5p . Same parameters as in 
Fig.El 



of the random telegraphic-like form of the fluorescence 
signal [a, |26|, being the appearance of dark periods, in 
which the atom does not emit photons, associated with 
the excitation of state |2), and hence giving rise to a line 
broadening of the ground state |g). In our case, this dy- 
namics is due to the mechanical effect of light, and indeed 
Tb scales inversely proportional to the localization of the 
atoms about the trap center. 

Let us now analyze the corresponding features of the 
quantum jumps in the spectrum of resonance fluorescence 
of the light emitted along the transition \g) -O- |2). Fig- 
ure S] displays the spectrum derived by numerically eval- 
uating the various terms to second order in the Lamb- 
Dicke expansion. As there is no excitation of state |2) 
at lowest order, all contributions to the spectrum are at 
second order in t]. We observe two features (also in this 
case the elastic peak is not reported): the blue sideband 
of the elastic peak with one large pedestal. This charac- 
teristic is independent of the angle of observation. The 
other sideband, corresponding to the transitions heating 
the atom, is hardly visible. From the magnifications in 
the insets, one can see that its height is two orders of 
magnitude smaller than the height of the blue sideband. 

We first discuss the asymmetry between the two side- 
bands. The spectra are evaluated at steady state of laser 



cooling where detailed balance of the heating and cool- 
ing rates between the phononic populations is satisfied. 
These rates are given by the transition rates connected to 
the scattering processes of laser cooling, multiplied by the 
population of the vibrational levels. For a single atomic 
dipole involved in the laser cooling, the strengths of the 
red and blue sidebands integrated over all directions of 
emission, is essentially given by these heating and cooling 
rates, which are equal. For multiple atomic transitions 
participating in the cooling, the motional sidebands re- 
flect only the scattering processes belonging to the con- 
sidered dipole, and thus their strengths can differ. In 
the competition of the two cooling mechanisms, the side- 
band cooling of the standing wave laser coupling transi- 
tion \g) -f-J- |2) dominates for the parameters used here. 
As a result, in the spectrum of this transition we observe 
a strongly pronounced blue sideband and an almost com- 
pletely suppressed red sideband. This asymmetry is bal- 
anced out by the other situation in the angle-integrated 
spectrum of the Doppler cooled transition. Our intuitive 
explanation is also supported by an analytical calcula- 
tion, which can be done at leading order in the satura- 
tion parameter s, and is reported in Appendix ID] The 
situation encountered here is reminiscent to the case of 
cavity cooling of a trapped atom [20| , where also two de- 
cay channels - the spontaneous emission of the atom and 
the cavity loss - exist. 

The pedestal appearing at the location of the blue side- 
band is an inelastic contribution to the spectrum, and 
arises from the broadening of the blue sideband transi- 
tion due to the coupling with state |1). Its central fre- 
quency is determined by the choice of the parameters of 
the laser coupling to state |2). An analytical calculation 
shows that it has the form 



^.M 



Td 



71/2 



w~WL2+i/ Ts 7f /4 + (cjl2 + i^ - w)^ 



(41) 



where 7^ = 571 /2 and which has been derived assuming 
s ^ 1. In the regime discussed in Sec. IIII Al in which 
quantum jumps are observable, it is hence dominant. In- 
terestingly, its height is proportional to the temperature 
of the atom, as one can check using the explicit form 
of the time scales To and Tg, Eqs. (|26t and (|25|) . and 
is also proportional to the degree of localization of the 
atom around the trap center, Ax^ ~ £,'^{n)- 

Wc finally discuss the dependence of the spectral sig- 
nals on the angles (j)i and (/)2 of the incident lasers and on 
the detector position V'- The ■(/'-dependence enters both 
in the dipole pattern of the emitted radiation, which is 
an overall factor, and in the individual spectral compo- 
nents. We focus on the latter ones. In the limit con- 
sidered here, the motional sidebands at w ~ ojli ± i^ are 
the only signals depending on ip. Their angle dependence 
is discussed in Ref. [22| and emerges from the quantum 
interference between photon scattering processes leading 
to an exchange of a vibrational excitation. For (j)i = 7r/2, 
i.e. when the driving laser is perpendicular to the axis 
of motion, the motional sidebands are only due to dif- 
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FIG. 4: Spectrum of resonance fluorescence of the light emit- 
ted by state |2). The sidebands are magnified in the insets. 
Same parameters as in Fig. (2] 



fusive processes, resulting in a cos '0-depcndence. The 
strength of the central peak in the spectrum of transition 
\g) -H- |1) is proportional to Tg^ ex cos^ 02, hence it de- 
pends on the orientation of the second laser driving the 
other transition \g) o |2). Similarly, our results show 
that 5'^^''(aj) ex cos^02: AH the spectral components of 
this transition are due to the mechanical effects of the 
second laser, and disappear for 02 = 7r/2. Moreover, 
the laser angles determine the mean vibrational excita- 
tion number (n), Eq. (jB7p . This quantity additionally 
scales the components of the spectra which are related to 
mechanical effects, as shown in App.ID] The spectral fea- 
tures connected with the quantum jumps are optimally 
observed when the wave vector of the standing wave is 
aligned with the atomic motion (02 = 0) and the de- 
tector angle "0 is at the maximum of the dipole emission 
pattern. In the regime considered here, the orientation of 
the running wave laser is less essential, best visibility we 
found for illumination along the axis of motion (0i — 0). 



IV. DISCUSSION AND CONCLUSIONS 

We have theoretically studied the resonance fluores- 
cence from a confined atom in a setup, where one could 
observe quantum jumps which are solely due to the me- 
chanical effects of light. In the limit in which the lasers 
driving the atom also cool the motion, we have derived 
the time scales of the bright and dark periods of the emit- 
ted intensity and the corresponding features in the spec- 
trum of resonance fluorescence. The latter is a pedestal 
appearing in correspondence of the elastic peak for the ra- 
diation emitted along the bright transition, whose height 
is proportional to the temperature and degree of local- 
ization of the atom about the trap center. 

We now discuss the experimental feasibility of this 
scheme. The setup we have considered assumes a V- 
configuration of levels, where one of the excited state is 
metastable. The atom is Doppler cooled along the op- 
tically allowed transition, while the transition which is 
metastable is driven by a standing-wave laser. In this 



work we assumed that the linewidth of this transition is 
much narrower than the trap frequency, so that the cor- 
responding laser essentially drives the red sideband. This 
is however not an essential requirement. The important 
condition is that the linewidth of the metastable state 
is significantly smaller than the rate of photon scatter- 
ing along the allowed, optical transition. A situation like 
the one here considered could be encountered in '*°Ca"'" 
ion, where the F-transition could be constituted by dipo- 
lar transition S1/2 — ^ ^1/2 , where the atom is Doppler 
cooled, and quadrupole transition 5*1/2 ~^ -D5/2, which is 
driven by a standing wave field, like the one of an optical 
cavity, as realized in (45[ . Since state 1^5/2 has a lifetime 
of the order of a second, this could be quenched by cou- 
pling to -P3/2, thereby not affecting the quantum jump 
scheme. 

While the motivation of this work is mainly to inves- 
tigate some fundamental properties of quantum jumps 
in resonance fluorescence, one could consider which out- 
looks on quantum engineering this work may offer. Be- 
yond the possibility to extract temperature and localiza- 
tion of an atom, one could consider a setup in which the 
quantum jump corresponds to the transition to a nonclas- 
sical state of the motion. In this case, this would allow 
one to post-selecting the quantum state, in the spirit of 
the proposal in Ref. [T2|. IJtI, lisj . 
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Appendix A: Eigensystem of d 

In this appendix we classify the left and right side 
eigenelements pj ^ and pj ' of the Liouville operator for 
the internal motion Ci, Eq. (|13p . obeying 

Cip^' = Xip^\ p^'Ci = Xip^\ (Al) 



A part of the eigensystem of £1 can be most easily cal- 
culated using the right and left eigenstates 

l±> = (^±|1> + ^l.9>) /^u;l + ny4 (A2) 

(^ = ("..±(11 + ^{g\\ l^ul + nl/A (A3) 



and 1 2) of the effective Hamiltonian 

H^P = H,^, + y(0) -^hY, |b-)0-| (A4) 



with eigenfrequencies 



w± 



and u)2 



^yJ±2Vi'-^T 



r!2 



(A5) 



, - i^. By using |2)(±| and |2)(±| together 
with their Hcrniitian conjugates as an ansatz, one finds 

p^ = |2)^, /5j^^-|±)(2| to Ai,± = (c.± - ^)A 

(A6) 

p^=^(2|, pI'^\2){±\ to A2,± = (c.2-^i)A. 

(A7) 

Apart from these four elements which describe the tran- 
sition between |2) and the dressed \g) -^ |1) transition, 
we further have 

pI- = |2)(2|, pI- ^ |2)(2| - ,^J_^(|,|2|i)(i| 



[T + kP]|5)(5|-mil)(sl+H.c.} 



(A8) 



to A. 



-72 with T 



71-272 ^"^ "> ~ 25i+i(7i-272) 



describing the decay of population in |2). The remaining 
four eigenelements are those of the two-level subsystem 
\g) O |1) with Rabi frequency J7i, detuning 5i and decay 
rate 71. These are 



Pi' = 1, Pi' = Pst to Ast = 
belonging to the internal steady state 

Pst^j^ {nl\i){i\ + iMst-nl)\g){g\ 



(A9) 



where 



{ni[2Si+iji]\l){g\+Ti.c.} 



Ust ^ if + 46l + 2nl, 



(AlO) 



(All) 



and the three eigenvalues AipLg, A^^g. From the latter 
ones, AlJ^Lg is real, the other two form a complex conju- 
gate pair at saturative driving. They represent the pos- 
sible transitions between the states of the dressed two- 
level subsystem and correspond to its MoUow-triplet. For 
(5i = simple analytic expressions for these eigenvalues 
and eigenelements can be found fig • 



Appendix B: Cooling 

The atomic motion is cooled by the mechanical effects 
of the two lasers. The motional state at the final stage of 



the coo ling dynamics can be estimated in perturbation 
theory [40, \^. At lowest relevant order in the Lamb- 
Dicke parameter, the dynamics of the population p„ = 
Tt {\n){n\g} in the vibrational state \n) is described by 
the rate equation 

Pn ^{n + l)(^l- + A2-)pn+l + n{Ai+ + A2+)pn-l 

(Bl) 

- [{n + l)(Ai+ + A2+) + n(Ai_ + A2-)] Pn (B2) 

where (n + l)(yli_|_ + A2+) and n(Ai_ + A2-) are the 
rates for population transfer between the levels n and the 
levels n + 1 and n — 1 respectively. The two terms Ai± 
and A2± entering these rates account for the mechanical 
effects of the lasers coupled to the transitions \g) <-^ |1) 
and \g) o |2). They are given by 



Ai± = 2Re{si{Ti^) + D} 
A2± = 2Re{s2(Ti^)} 



(B3) 



where 



S]{^) 






dre'^^Tr 



|V^0-),£..^0)^^^| 



71 



D = vtm^Tr{\l){l\p,,} 



(B4) 



with 






(1) 



(2) 



= ih — k-[ cosd 
2 

= —h—-k2 cos 



■i[|l)(5|- 
02[|2)(.g| 



l.9)(l|] 
H5>(2|] 



(B5) 



and Pst being the internal steady state from Eq. (jAlOp . 
The motional state at the final stage of cooling 



1 



Pst 



1 



1 



(B6) 



follows from the rate equation (jB2|) for p„ == 0. The av- 
erage number of vibrational excitations can be expressed 
in the form 



niWi + n2W2 
W 



(B7) 



where Uj = Aj^/{Aj^ — Aj^) with j = 1, 2 is the excita- 
tion number corresponding to the situation in which only 
the transition \g) -?->■ \j) participates in the cooling pro- 
cess, and Wj — Aj_— Aj^ are the corresponding cooling 
rates. The total cooling rate which takes into account 
both transitions is given by 



W = Wi+W2. 



(B8) 



The steady state of the motional dynamics results from 
the competition of the mechanical effects of the photons 
scattered along the two transitions. 

Simple expressions for rij and Wj can be found when 
the parameters are set for Doppler cooling on the tran- 
sition |1) o \g) and sideband cooling on the transition 
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FIG. 5: Mean phonon number (n) for the different cooling 
schemes as function of 62- The dashed hne shows the result- 
ing phonon number ni of the Doppler cooling only, while the 
dotted line corresponds to the mean phonon number if there 
was only sideband cooling at the standing-wave transition. 
The solid line shows the mean phonon number from Eq. (IB7[l 
resulting from the competition between the two cooling mech- 
anisms. Same parameters as in Fig. [2] where 82 was chosen 
to optimize the phononic excitation number (n). 



|2) -n- \g), i.e. when 71 > v, f2i, 5I2 > 72, Si = 71/2 and 
62 — Rew_ -I- v. The terms tt-i and Wi are equal to the 
standard results of laser cooling of a trapped two level 
atom and are given by 



Til 






Wi w 2r]l- 



1 + -^ 
cos-^ 01 

M'^ cos^ (j)i 



(B9) 



The cooling on the transition |2) o \g) is influenced by 
the dynamics on the other transition which is relevant at 
zero order, see App. [Xj In particular, the corresponding 
terms n2 and W2^ are found by taking into account the 
light shift Rew„ ~ ~j-^ of the ground state due to the 
coupling to the first laser. Together with the light shift 
the ground state acquire a finite line width 7_ 
the corresponding results can be expressed as 



tH- and 

271 



n2 
W2 



(72+7-) 
16t/2 

772^^2 COS^ 02 

72+7- 



(BIO) 



Sm("^lA*st|"T-)Pm(i) in the basis of the vibrational eigen- 
states |?7i) and expand 



^™W = ^[. 



-iO„ 



(X2,n' |X2n) (X2n |to) \g) (g| (m|x2,«') 



CT.cr'— ± 



''""' (x<T',n" |XrT,«) (x<T« |to) \g) {g\ {m\xa 



(CI) 



in terms of the right eigenstates |x±,n), 1x2, «) and left 
eigenstates (x±.n|, (X2,n| of the effective Hamiltonian 
ffcff, Eq. ([^ . with eigenvalues huj±n, tujj2n- We also 
defined the slowly decaying rates ^nn' = '^2n — ^2n' ^^'^ 
the fast decaying rates ujg-a'nn' = ^an — '^a'n' ^^ abbre- 
viations. In the last step we have neglected the mixed 
terms which include rapid and slowly decaying frequen- 
cies: They are only small corrections to the second line 
in Eq. (ICTj) . 

To proceed, we evaluate Eq. (ICip pcrturbatively in the 
Lamb-Dicke parameters rjj, Eq. (jO]). To this end, we ex- 
pand the effective Hamiltonian, 



r(off) 
^0 



H 



(off) 



its eigenfrequencies, 

w±„ = 

W2n = Wn 

and the right and left eigenvectors 

|X2,n) = |X2,n)o+ |X2,n)l 



,,(±") 
(2,n) 



,(±") 
(2,n) 



(X2,«| =o(X2,«| +l(X2,«| + 
IX±,n) = ••■ 



(C2) 

(C3) 
(C4) 

(C5) 
(C6) 
(C7) 



(off) 



in increasing order in r/j. In lowest order, Hq is a sum 

of the external part Hcxt and the internal part i?j°^ 
from Eq. (|A4p whose right and left eigenstates |±), |2) 
and (±1, (2 1 with their eigenvalues u!± and 102 are listed 
in App. El Eqs. (|X3])-(|M1)- Moreover, we have 

|X2,«)o = |2)|n), \x±.n)o ^ \±)\n) (C8) 

and analogous expressions for the left eigenstates. The 
higher order corrections we find using standard tech- 
niques of perturbation theory. We now return to Eq. (|Cip 
and evaluate 



Pm{t) 



A^' 



-''''\i{X2,n'\9n)\^ 



T.-"'^ 



-<')* 



a\a'){a\g){g\a'l). (C9) 



Appendix C: Evaluation of the bright and dark time 

scales 



In order to find expressions for the bright and 
dark periods, we start from Eq. (PT|) . rewrite P{t) = 



up to second order in the Lamb-Dicke parameter. It turns 
out that we only need to calculate the first order correc- 
tion 



i(X2,n| =?y2cos02 (v^(n - 1| + VriTT(n+ 1|) 

hn2/2 



x(g|- 



huj, 



(2,n) 



H, 



(off)- 



(CIO) 
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Using expressions ([£9]) and (|C10|) in Eqs. ([23]) and ([24]) 
finally yields the expressions Tg and T^ from Eqs. (^5)) 
and (f26| when taking into account the separation of the 
time scales |Ima;±| 3> - S> |Ima;2|. In Eq. (|25p we use 
the abbreviations 



T(o) = -z ^ 



{<j\a'){a\9){g\a') 



1 



fT.c' — lb 



(w^ 



7i(1|pS*|1) 



£=±1 



i{n)+Si,i) 



E 

(T = ± 



(skX'Tlff) 



W2 



£iy 



(Cll) 



(C12) 



The time T^^' is the average waiting time between two 
photon emissions of the atom at zero order in the Lamb- 
Dicke expansion and F determines the temperature de- 
pendence of the bright period. An expansion in the sat- 
uration parameter s of this quantity gives 

16(n) f_ . ^472^ 
S7i. 



T- 



ih^ 



1 



(C13) 



Appendix D: Evaluation of the spectrum of 
resonance fluorescence 



In this appendix we sketch the derivation of the analyt- 
ical expressions for several components of the spectrum 
of resonance fluorescence. The central peak of transition 
\g) -i-> |1) can be calculated starting from Eq. (15^ with 
A^ = —72 and Ae = 0. Wc find 

Sg^ (u;) = Re-^— l^^^a) (Dl) 



z(w-a;Li)-A^ "^ 



where only 



F(i) = Tv{Dlv^Dofy^n + '^HDlv>:Dopf} 



A=A-, 



(D2) 



of the ten possible second order combinations of the 
internal traces -F2(A), Eq. ([ST]) , contribute. The re- 
quired techniques for evaluating the remaining traces can 
be found, e.g. in Refs. [19|, l23j. It turns out, that 
Tt{DIV2DoPq} cx 7^^ dominates. Constraining our- 
selves on this term and exploiting 72 ^71, we find 



p(i) 

cp 



2 „„„2 i, ^'2 



272 



(D3) 



where 



Pi{n))^ReY,{g\a){a\pst\9) 



{2{n) + l){j2 + Xia) - ii^ 



(72 



Al.)2 



_ (D4) 

Here, \a) and \a) denote the left and right eigenstates of 
the effective Hamiltonian, (|A4|) . and \j± are eigenvalues 
of the Liouville operator £1 discussed in App. [Xj Taking 
into account the conditions for Doppler cooling and side- 
band cooling at the corresponding atomic transition, we 
obtain for the peak's weight function 



p(i) 

cp 



772 cos^ (f>2 (n) 



7172 



1 -2s 



272 
s 71 



(D5) 



up to first order of saturation parameter s. 

We now turn to transition \g) o |2). Here, three spec- 
tral signals of second order are observable: Two motional 
sidebands of the elastic peak. 



si^H^Re 



1 



A^), 



i(cl; -CJL2) ±«i' + 7sb 



and an inelastic contribution, 



:(2) 



^i\:» = Re 



i{uj - LJL2) + Ai- 



v(2) 

inel' 



(D6) 



(D7) 



belonging to the eigenvalue Ai_ of £1. In Eq. (|D6I) we 
introduced the width of the motional sidebands, 7sb = 
W/2, given by the real part of the second order correction 
A2 and corresponding to the cooling rate W, Eq. (|B8I) . 
of the laser cooling process [42 1. 

In order to characterize the form of the spectral sig- 

nals, we need to evaluate the weighting terms F^^'{±ii') 

(2) 
and -Find. ^'^^ ^^^ motional sidebands, only the term 

Ti{DlViDQpf} contributes and yields 



^sb («M='72COS2 02 



nU{n)+5e,-i) 



E' 



{g\cr){a\pst\g) 



Aicr — Ae 



(D8) 



The real expression Eq. (|D8[) indicates that the motional 
sidebands are pure Lorentzians and independent from the 
detector position angle ip. An expansion in the saturation 
parameter s under the same conditions as for Eq. (jDSp . 
reveals 






16l^2 



(1 - s) (D9) 



(•21. N 9 9 , 4r2o(n) / 472 

^i' M='?2COS2 02^V 



1 



7fs- \ S7i^ 
It remains to address the term of the inelastic peak. 



(DIO) 



F; 



(2) 



2 2 I '^'2 



incl 



-r;^cos"<^2:P(ffh)'x 
272 

{9\<y){^\Pst\9 y ^ \2 , 1 (Dii) 



a=± 



AL + -2 



from Eq. (|D7I) . To arrive at Eq. (JDIII) , we have taken into 

account only the dominating term TrjIJjT'o-DoPl'} °^ 
7^ . Analogously to the previous treatments, an expan- 
sion in s yields 



i;.(2) 2 2 , 2f7^ (n) 

^i„cl=%COS 02-— — 

7172 s 



1_ lll^it 

S7i *2 



(D12) 



When using the approximations for the dark and bright 
period To and Tg from Sec. IIII Al we eventually find the 
expressions PO]) and (|4ip . in the case 72 <C S71/2. 
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